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FATIGUE  CRACK  TIP  STRESS  AND  STRAIN  RANGE  DISTRIBUTIONS,  THE  ENERGY 
FOR  CRACK  PROPAGATION  AND  THE  EFFECT  OF  WATER  VAPOR  ON  THESE  FACTORS 

D.  L.  Davidson  and  J.  Lankford 

Introduction 

A knowledge  of  the  stress,  strain,  and  energy  distributions  around 
the  tip  of  a fatigue  crack  would  be  helpful  In  understanding  the  response 
of  the  material  to  the  remotely  applied  loading  conditions,  both  for 
the  condition  of  environmental  interaction  and  without.  Then  from  the 
material  response,  it  might  be  possible  to  either  predict  the  rate  of 
fatigue  crack  propagation  or  improvise  a method  of  retarding  it. 

To  this  end,  Davidson  and  Lankford have  utilized  the  propensity  of 
low  carbon  steel  to  form  a dislocation  cell  structure  near  the  crack  to 
calculate  the  energy  required  to  propagate  the  crack.  From  this  calculation, 
the  spacial  distribution  of  absorbed  energy  around  the  crack  may  also  be 
determined.  It  is  the  purpose  of  this  paper  to  extend  the  use  of  the 
experimentally  measured  dislocation  cell  structure  distribution  to  estimate 
the  stress  distribution  in  the  near  crack  region. 

Criteria  for  dislocation  cell  formation 

Dislocation  cellular  structures,  or  subgrains,  form  within 
the  grain  structure  of  low  carbon  steel  upon  the  attainment  of  a critical 
dislocation  density.  The  number  of  cycles  N required  for  subgrain  formation 
is  dependent  on  the  level  of  shear  strain  per  cycle.  The  details  of 
subgrain  formation  are  different  on  the  surface  and  in  the  interior,  but 
once  the  critical  dislocation  density  is  achieved,  and  subgrains  form, 
continued  cycling  causes  the  subgrain  size  to  decrease  to  an  approximately 
equilibrium  value. 

For  a smooth  specimen  cycled  between  controlled  limits  of  plastic 
strain,  giving  a plastic  strain  range  Aep,  the  resulting  stress  range, 

Aa,  at  first  changes  relatively  fast  and  finally  reaches  an  approximately 
equilibrium  value  (but  oscillations  of  unknown  origin  sometimes  do  occur). 
Attainment  of  an  equilibrium  stress  for  this  loading  condition  is  thought 
to  coincide  with  the  formation  of  an  equilibrium  subgrain  size. (2,3) 

For  a smooth  specimen  cycled  between  controlled  limits  of  stress, 
resulting  in  a range  of  applied  stress,  Ao,  the  plastic  strain  range  at 
first  increases  rapidly  and  finally  reaches  an  approximate  equilibrium 
value,  although  subsequent  small  oscillations  in  this  value  have  been  found(^). 
Attainment  of  an  equilibrium  ACp  is  thought  to  coincide  with  the  formation 
of  an  equilibrium  subgrain  size. 

Likewise,  the  work  per  cycle,  as  measured  by  the  area  within  a plot 
of  Ao  vs  A£p,  increases  rapidly  with  cycling,  and  finally  attains  an 
approximately  equilibrium  value  which  is  believed  to  coincide  with  an 
equilibrium  subgrain  size.  It  is  this  last  effect  which  allowed  calculation 
of  the  energy  of  crack  propagation  to  be  determined  from  the  measured 
subgrain  size  distribution  around  a fatigue  crack. (^)  Subgrain  sizes  may 


r ■ ■ ' 

II 

, i 

I also  be  used  as  a measure  of  stress  distribution  due  to  the  attainment  of 

; an  equilibrium  value  of  stress  under  strain  controlled  cycling,  and  may  be 

used  to  measure  strain  distribution  if  a way  of  summing  AEp  over  N can 
, be  found  to  obtain  total  or  accumulated  strain. 

Crack  Tip  Subgrain  Distribution 

The  distribution  of  subgrains  in  the  wake  of  a fatigue  crack  has 
been  measured  for  a range  of  cyclic  stress  intensities  (AK) , both  for 
crack  propagation  in  an  inert  environment  (nitrogen  with  1 ppmv  water 
vapor)  and  in  an  aggressive  environment  (12,000  ppmv  water  vapor  or  50% 
R.H.).  The  subgrain  distribution  was  determined,  using  the  electron 
channeling  contrast  technique  (5) , and  resulted  in  the  finding  of  a 
linear  relation  between  subgrain  size,  d,  and  distance  from  the  crack, 
r.(6) 

d = A + Br  (1) 


There  is  statistical  scatter  in  the  measurement  of  A and  B from 
specimen  to  specimen.  The  reasons  for  this  scatter  are  partially  known. 
Clearly,  there  is  an  interaction  between  the  crack  tip  stress  field  and 
the  local  grain  orientations  near  the  crack  which  has  some  random  nature. 
Microscopically,  the  crack  deviates  from  growth  perpendicular  to  the 
loading  axis,  thereby  causing  a changing  mix  of  mode  I and  mode  II 
opening  as  the  crack  advances,  which  is  known  to  alter  the  crack  tip  stress 
field.  The  data  currently  available.  Figure  1,  for  the  intercept  A, 
slope  B,  and  extent  of  subcell  formation  r^,  suggests  that  the  presence  of 
water  vapor  causes  the  subgrain  size  near  the  crack  plane  to  undergo 
alteration  as  AK  increases,  whereas  an  inert  environment  does  not. 

For  both  environments,  increasing  AK  causes  the  slope  B to  decrease, 
with  B for  the  water  vapor  environment  decreasing  faster  than  for  the  dry 
environment.  More  data  is  required  for  a refinement  of  quantitative 
details  of  these  observations,  but  the  trends  are  clear. 

Equation  1 also  correlates  subgrain  size  with  distance  ahead  of  the 
crack  tip,  but  data  in  this  region  is  the  result  of  a decreased  number  of 
cycles,  and  may  be  the  result  of  a nonequilibrium  subgrain  size.  The 
change  in  subgrain  size  also  occurs  over  a much  shorter  distance  in  the 
crack  tip  region  making  the  slope  difficult  to  determine,  except  by 
averaging  (i.e.,  between  the  subgrain  size  at  the  crack  tip  and  the  grain 
size) . The  size  and  shape  of  the  subgrain  forming  region  near  the  crack 
tip  is  also  known,  as  illustrated  in  Figure  2. 

Correlations  of  Cyclic  Stress  Range  and  Subgrain  Size 

The  number  of  published  research  results  relating  subgrain  size  to 
stress,  strain,  or  energy  dissipation  is  not  large,  but  several  pieces 
of  excellent  work  have  been  accomplished.  A summary  of  the  available 
results  is  shown  in  Figure  3.  The  work  of  Lawrence  and  Jones(^),  by  far 
the  most  complete,  is  complemented  by  that  of  Chopra  and  Gowdaw),  and 
the  one  experimental  point  from  Davidson  and  Lankford. It  is  interesting 
that  the  results  of  Herz,  Mughrabl  and  Wilkens,^®^  while  not  quantitatively 
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Figure  1.  The  Effect  of  Cyclic  Stress  Intensity  Factor  on  Various 
Experimentally  Derived  Variables 
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Figure  3.  Cyclic  Stress  Range  Vs  Subgrain  Size.  Also  Shown  Is 
Data  From  Monotonlc  Loading. 


equivalent,  have  about  the  same  slope.* 


The  results  for  tensile  testing  have  neither  the  same  slope,  nor  the 
same  range  of  subcell  size,  except  for  the  work  of  Keh  and  Weis sman, which 
is  in  the  same  range  as  part  of  the  results  of  Chopra  and  Gowda. 

The  Subgrain  size  - dislocation  density  relation 

Questions  of  interpretation  arose  during  the  examination  of  data 
for  inclusion  in  Figure  3.  By  examining  the  relation  between  dislocation 
density  and  subgrain  size,  some  of  these  questions  can  be  resolved. 

Experimentally,  the  relation  between  subgrain  size  and  dislocation 
density  has  been  found  by  many  investigators  to  be 


where  K©  is  a constant  (probably  material  dependent).  From  theoretical 
considerations,  dislocations  are  expected  to  rearrange  themselves  into 
periodic  clusters  once  a critical  density  has  been  achieved.  The  mathematical 
model  for  this  process  was  developed  by  Holt, (5)  and  it  predicts  the  rela- 
tion of  equation  (2) . This  relation  will  thus  be  considered  as  a basic  physi- 
cal relation. 

For  tensile  deformation,  a relation  may  be  derived  between  the  shear 
strain  and  the  dislocation  density^^O^ 


Y - vT 


(3) 


where  = a material  constant,  and  Yq  is  a constant  of  integration,  which 
physically  may  be  zero.  By  combining  equations  (2)  and  (3)  the  relation 


Yd  = K2 (constant) 


(4) 


is  obtained. 

(2) 

Abdel  - Raouf,  et  al  found,  through  strain  rate  change  experiments, 

that 


d = 


(5) 


which,  by  rearranging,  gives 


* The  results  plotted  are  saturation  stress  Og,  even  though  Herz,  et  al,  find  the 
athermal  component  of  the  cyclic  flow  stress  Ogg  to  correlate  better.  For 
the  strain  rates  and  strain  amplitudes  relevant  to  crack  propagation,  Ogg 
is  a large  portion  of  Og,  and  should  correlate  with  subgrain  size  nearly  as  well. 


Equations  (4)  and  (6)  predict  different  relations  between  strain  and 
subgrain  size,  which  may  be  due  to  monotonic  vs  cyclic  loading  or  to  the 
experimental  technique  used  by  Abdel-Raouf,  et  al  in  their  data  collection. 

By  using  the  Seeger  theory  of  the  dependence  of  stress  on  temperature 
and  strain  rate  for  cyclic  deformation,  as  expounded  by  Mughrabi,  et  al,'^^^ 


- ■'SG  ' 


where  Tgg  is  the  athermal  component  of  the  stress,  and  Xg  is  the  saturation 
stress,  and  substituting  in  equation  (5),  a cyclic  stress-strain  curve  of 
the  form 


s SG  3 


(^) 


is  obtained,  which  fits  well  with  the  cyclic  stress-strain  curve  of 
Landgraf  (12) 


Ao  = 2K' 


where  (experimentally)  n'  = 0.27  is  the  cyclic  work  hardening  coefficient, 
and  K'  = 220  for  this  Fe  - .05C. 


Herz,  et  al  ' find  that  there  is  a correlation  between  dislocation 
density  and  subgrain  size  and  the  athermal  cyclic  stress  Xgg  for  pure 
iron  single  crystals.  They  find 


1.6  , , -1-4 

p « and  d « Tgg 


which  leads  to  the  relation 


dp  “ constant 


This  relationship  does  not  fit  equation  (2) . By  reinterpreting  the 

data  of  Herz,  et  al(8)  so  that  both  p and  d are  proportional  to  the  saturation 

stress  Tg,  then  the  following  relations  are  obtained 


d - 2 X 10^  T 

s 


J 


which  compares  well  with  equation  (12) 

The  conclusions  drawn  from  the  above  analysis  are: 

a)  The  functional  form  of  equation  (13)  or  (14)  is  correct. 

(3) 

b)  Our  analysis  of  the  data  of  Herz,  et  al  is  consistent  with 
both  the  data  of  Figure  3 and  the  physical  relation  given  by 
equation  (2). 

c)  There  is  an  inconsistency  in  the  correlation  given  by  Abdel-Raouf, 

et  equation  (5). 

The  self-consistency  of  the  above  argument  fits  with  the  experimental 
correlations  made  previously  between  work  per  cycle  and  subcell  size  (Ref  1, 


8 


Appendix  C) . It  also  indicates  a consistency  between  the  mechanism  of 
subgrain  formation  in  tensile  and  cyclic  deformation,  implies  that  the  concept 
of  Similitude (10)  is  correct,  and  reaffirms  the  validity  of  equation  (2). 

The  stress  range  - subgrain  size  correlation 

The  results  we  have  chosen  to  use  for  the  calculations  made  in 
this  paper.  Figure  3,  are  shown  as  the  heavy  black  line  connecting  the 
yield  stress-grain  size  point  for  our  material,  and  passing  through  the  data 
of  Chopra  and  Gowda,  part  of  that  of  Lawrence  and  Jones,  and  the  one 
experimental  point  obtained  for  our  material.  In  effect,  we  have  normalized 
the  data  of  Lawrence  and  Jones  and  Herz  and  Mughrabl  to  our  grain  size  material. 
The  reason  for  using  the  grain  size-yield  stress  point  is  that  subgrains  are 
known  to  begin  forming  at  the  yield  stress.  This  is  illustrated  in  Figure 
4,  which  shows  the  material  on  either  side  of  a Luders  front.  As  may  be 
seen,  the  passage  of  the  Luders  front  results  in  subgrain  formation. 

Strain-rate  effects 

Mughrabi,  Herz  and  Stark^^^^  have  clearly  shoim  that  strain  rate 
is  an  important  variable  in  assessing  the  deformation  properties  of  iron 
single  crystals.  Similar  work  on  rate  effects  for  poiycrystalline  iron 
has  been  done  by  Abdel-Raouf,  Behnam  and  Plumbtree'^’^) , which  shows  the 
importance  of  this  factor  although  their  results  are  inconsistent  with 
those  of  Herz,  et  al.  A comparison  of  the  present  experimental  results, 
using  a loading  rate  of  5 Hz  (a  crack  tip  strain  rate  of  approximately 
10  sec“^)**  with  those  previously  reported  for  the  same  Fe-.05C  steel, (^®) 
for  which  the  loading  rate  was  152  Hz**  (a  crack  tip  strain  rate  of  approxi- 
mately 3 X 102  sec~l)**.  Figure  5,  reveals  a difference  in  the  extent  of 
plasticity  along  the  crack  wake  and  a difference  in  the  effect  of  AK  on 
this  plasticity. 

Because  of  these  conflicting  data  on  rate  effects,  the  data  in 
Figure  3 have  not  been  corrected  to  one  value  of  strain  rate.  For  the 
strains  near  the  crack  tip  (~1) , saturation  stress  range  is  not  thought 
to  be  greatly  affected  by  a rate  change  of  10;  thus,  the  correlation  shown 
in  Figure  3 is  not  expected  to  be  much  different  between  5 and  50  Hz, 
which  covers  most  of  the  data  used  in  the  correlation. 

Crack  Tip  Stress  Range  Distribution 

Combining  the  data  in  Figures  1 and  3,  through  equations  1 and  14, 
the  stress  range  distribution  near  the  fatigue  crack  may  be  derived.  Stress 
range  has  been  computed  around  the  crack  tip  for  the  dry  nitrogen  environment 
case  for  two  values  of  AK,  and  for  both  the  dry  nitrogen  and  wet  air  cases 
for  two  values  of  AK.  The  equation  for  stress  range  is 


**  strain  rate(ll)  = 2 x frequency  x strain  range.  Crack  tip  strain  range  was 
assumed  to  be  1. 


Luders  Front 


Figure  4.  Formation  of  Subgrains  in  Fe  - .05C  Steel  Due  To 
the  Passage  of  a Luders  Front  Upon  Yielding.  The  Front 
Moved  From  Left  to  Right,  With  Material  On  the  Left 
Having  Formed  Subgrains. 
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Figure  5.  Correlation  Between  Plastic  Zone  Size  and  Cyclic  Stress 
Intensity  Factor  For  Cracks  Propagated  In  Humidity  Air 
(Approximately  50%  RH)  Showing  the  Effect  of  Frequency 

of  Loading. 
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or  Che  distance  from  the  crack  tip  may  be  normalized  by  the  factor  propor- 
tional Co  the  plastic  zone  size  (AK/Oy)^,  giving 


r - ^ (18) 

m/a^r 

80  that 

An  - “-262 

— - 2.26(A  + Br)  (19) 
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The  effect  of  changing  AK 

The  stress  range  distribution  around  the  crack  tip  has  been 
computed  for  AK  » 8 and  12  MN/m^/2.  71,0  functions  used  to  calculate  the 

distribution  are  shown  In  Figure  6.  These  were  derived  from  the  observation 
Chat  plastic  zone  shape  Is  described  approximately  by  Che  plane  strain  clastic 
zone  calculation  as  modified  by  experimental  observation  near  0 0'’.(^) 

The  AK  dependence  of  the  plastic  zone  boundary  dimension  rj^  was  measured 
experimentally  as  AK^*^  at  6 - 90°,  and  the  values  at  the  other  angles  have 
been  assumed  to  have  the  same  AK  dependence. 

Using  the  above  observations  and  assumptions,  the  yield  stress 
normalized  stress  range  distribution  Is  found  to  be  as  Is  shown  In  Figures 
7 and  8.  The  distribution  of  stresses  Is  a smooth  function  having  the  form 
given  by  equations  (17)  or  (19)  which  approaches 


Ao 
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2.26 


(b7)-262 


(20) 


in  the  outer  part  of  the  subcell  forming  region. 

The  curves  cannot  quite  be  normalized  by  equation  (18)  because  Che 
data  of  Figure  1(c)  have  a AK^*^  dependency  Instead  of  Ak2.  thus  the  curves 
of  Figure  7 cannot  be  made  to  lie  one  on  the  other  using  T as  the  spaclal  variable. 

I 

The  effect  of  changing  environment 

Since  Che  subcell  size  near  the  crack  plane  Is  decreased  at  low 
AK  for  the  water  vapor  environment,  Che  Implied  cyclic  stress  range  is 
also  decreased,  and  the  stress  range  distribution  Is  changed,  reflecting 
the  fact  that  the  extent  of  plasticity  Is  related  to  AK  as  rj^  «Ak2*5 
(0  >•  90°,  AK  ■ up  to  13  MN/m3/2).  The  distribution  of  stress  range  around 
the  crack  tip  Is  shown  In  Figures  9 and  10,  In  comparison  to  Che  dry 
nitrogen  environment  for  the  same  AK.  At  the  higher  AK  (12  MN/m3/2), 
the  stress  range  distribution  becomes  nearly  equal  for  both  environments, 
and  appears  much  as  shown  In  Figure  7 . 
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Figure  6.  Angular  Relations  Between  Plastic  Zone  Size  and  Cyclic 
Stress  Intensity  Factor  At  the  Crack  Tip  For  Both  Dry 
Nitrogen  and  Wet  Air  Environments. 
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Figure  7.  Stress  Range  Distribution  Around  the  Fatigue  Crack  Tip 
For  a Crack  Propagated  in  Dry  Nitrogen  Showing  the  Effect 
of  Stress  Intensity  Factor. 
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Figure  9.  The  Effect  of  Environment  on  the  Crack  Tip  Stress  Range 
Distribution  At  AK  ■ 8 MN/m^/2 
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Figure  10.  The  Effect  of  Environment  On  the  Angular  Distribution  of 
the  Stress  Range  - Distance  Relation 
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The  rate  of  crack  propagation  has  very  accurately  been  measured  for 
the  water  vapor  and  dry  nitrogen  environments,  and  Is  shown  In  Figure  11. 
Only  the  resulting  statistically  fit  curves  are  shown  because  of  scatter 
In  the  data. 

Crack  Tip  Strain  Range 

Comparison  of  equation  (14)  with  that  of  the  cyclic  stress-strain 
curve,  yielded  the  relation  given  as  equation  (15)  , namely 


2 

A + Br 


(21) 


At  r » 0,  Aep  = 2 (dry  nitrogen) , Irrespective  of  AK,  and  Aep  - 0.53 
(wet  air,  AK  = 8 MN/m3/2),  jn  the  outer  portion  of  the  subcell  forming 
region,  Aep  “ ^ for  both  environments. 

The  Energy  Required  for  Fatigue  Crack  Propagation 

Subgrain  formation  by  passage  of  the  crack  may  also  be  used  to  cal- 
culate the  energy  required  to  create  a unit  area  of  new  fatigue  crack.  The 
methodology  for  this  calculation  has  been  fully  explained and  is  generally 
summarized  by  the  following  relations 


where 

W^(r)  » '^co*^  d = A + Br  (23) 

“ experimentally  determined  constant 
n * 1 

and 


N(r) 


dL 

da/dN 


(24) 


where  dL  Is  the  dimension  of  the  material  ahead  of  and  parallel  to  the 
crack  plane  in  which  the  subgrains  form. 

The  shape  of  dL(r)  has  been  approximated  as 

N(r)  - + Nj^(r)  o < r < r' 

N(r)  - N r'  < r < r^ 

where  r^^  Is  the  plastic  zone  size. 
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This  methodology  has  been  applied  to  the  data  thus  far  presented  for 
both  the  dry  nitrogen  and  wet  air  environments.  For  values  of  AK  above  20, 
the  slope  of  the  da/dN  vs  AK  Is  assumed  to  be  4.0  for  both  environments.  A 
"limiting  case"  calculation  may  be  made  for  AK  above  this  value.  To  make 
this  calculation  the  following  assumptions  have  been  made: 

<» 

2 

(2)  (plastic  zone  size)  « AK 

(3)  slope  B “ AK~^ 

The  data  of  Figures  1(c)  and  5 give  validity  to  assumption  (2), 
and  assumption  (3)  Is  reasonable  from  Figure  12. 

Using  these  assumptions,  the  "limiting  case"  gives  a value  of  fatigue 
crack  propagation  energy  Independent  of  AK.  For  values  of  AK  » 8 to  13  MN/m^/2,for 
the  dry  nitrogen  environment,  experimentally  derived  values  of  the  exponents 
and  of  B were  used.  For  AK  between  13  and  20,  Intermediate  or  "transition" 
values  were  used,  which  are  given  In  Table  I. 


TABLE  I 
Exponents 


AK 

da 

(MN/m3/2) 

dN 

’^l 

B 

8 

4.4 

1.6 

.22 

10 

4.4 

1.6 

.16 

12 

4.4 

1.6 

.12 

14 

4.3 

1.7 

.096 

16 

4.2 

1.8 

.071 

18 

4.1 

1.9 

.071 

For  both  environments, 
shown  In  Figure  13  and  Table 

the  resulting 
II. 

energy  values 

calculi 

TABLE  II 


Energy  for  Fatigue  Crack  Propagation 
(10^  Joules/meter^) 


Dry  Nitrogen 


Wet  Air 


Due  to  Environment 


8 

17.1 

1.5 

15.6 

! 

10 

12.5 

2.4 

10.1 

i 

1 

12 

9.9 

4.1 

5.8 

1 

14 

10.8 

9 (est) 

1.8 

16 

14.1 

14  (est) 

0 

18 

21.1 

LIMIT 

25 

Parametrically,  equation  (22)  has  the  form 


+ C2AK"'^ 


where  the  exponents  p and  q are  dependent  on  the  values  of  those  shown 
in  Table  I,  The  reason  Region  I (dry  nitrogen)  of  Figure  13  has  a negative 
slope  is  because  the  slope  of  da/dN  vs  AK  has  a slope  greater  than  4; 
conversely,  for  the  wet  air  environment  the  exponent  is  3.63  and  the  slope 
of  the  curve  in  Figure  13  is  strongly  positive.  The  magnitude  of  W.j.  in 
equation  (25)  is  also  dependent  on  Cj  and  C2  which  are  proportional  to  B; 
thus,  the  existence  of  Region  II,  or  transition  region  between  Region  I, 
the  threshold  region,  and  the  limiting  region  depends  very  strongly  on  the 
Interrelationships  between  the  actual  values  of  the  constants  chosen. 

The  previous  calculation  of  crack  propagation  energy^^^  found  for  dry 
nitrogen  Wf  = 17.3  x 10^  j/m^  and  for  wet  air  Wx  = 5.7  x 10^  J/m^,  both  at 
AK  ■ 10  MN/m^/^.  The  values  shown  in  Table  II  are  somewhat  different, 
and  reflect  greater  accuracy  in  the  determination  of  crack  propagation  rate. 

The  Effect  of  Hydrogen  on  Flow  Stress  and  Fracture 

The  effect  of  hydrogen  on  the  mechanical  properties  of  iron  and  iron 
alloys  has  been  the  subject  of  considerable  recent  effort.  The  work  of 
Karpenko,  Yarmchenko  and  Shved'^^'  showed  that  for  Armco  Iron  (7  x 10"^  wt  X C) , 
hydrogen  charging  Increased  the  yield  strength  from  296  to  336  MPa  when 
samples  were  thin  enough  to  obtain  hydrogen  charging  completely  through  the 
specimen  thickness.  Likewise,  the  tensile  ductility  was  reduced  by  hydrogen 
charging  during  testing.  Asano  and  Otsuka^^^^  examined  a decarburlzed 
mild  steel  (5  x 10“^  C)  and  found  essentially  the  same  results  as  Karpenko, 
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Figure  13. 
Cyclic 


Energy  Required  For  Crack  Propagation  As  a Function  of 
Stress  Intensity  Showing  the  Effect  of  Environment 


et  al,  at  least  qualitatively. 


Using  ultra  high  purity  iron  (3 
found  that  hydrogen  charging  resulted 
contrary  to  the  above  cited  results, 
work,  also  on  high  purity  iron,  which 
to  hydrogen. 

(18) 

Lunarska,  et  al^  ^ have  measured  by  torsion  pendulum  the  shear  modulus 
of  high  purity  iron  (5  ppm  nonmetalllc  impurities)  and  concluded  that  hydrogen 
charging  decreases  the  shear  modulus,  and  that  this  is  related  to  a reduction 
in  the  lattice  cohesive  strength. 

The  conclusions  to  be  drawn  from  these  studies  are  that  the  carbon 
content  of  the  iron  is  an  important  factor  in  whether  the  flow  stress  is 
Increased  or  decreased  by  hydrogen.  For  Fe  - .05C  steel,  then,  hydrogen 
charging  is  expected  to  increase  the  flow  stress.  Quantitatively  the 
effect  should  be  about  that  found  by  Karpenko,  et  al,  or  about  a 13%  Increase, 
assuming  saturation.  Hydrogen  charging  also  lowers  the  tensile  ductility  of 
steel,  probably  through  its  reduction  in  cohesive  strength. 

Interpretation  of  Results 

The  shapes  of  the  crack  tip  subcell  forming  regions  observed  for  this 
material  are  very  much  like  those  graphed  by  Schljve'^’^  who  derived  as}rmpotic 
solutions  for  the  ellptic  hole-in>«-sheet  solution  of  Tlmoshenka  and  Goodler. 
The  plain  strain  solutions  for  the  shear  stresses 

^ ^1  ~ ^3  j ^ _ ^2  ~ ^3 

^2  * 2 And  ” 2 


X 10  C),  Klmura,  Hatsul  and  Horlya^  ^ 

in  a decrease  in  flew  stress, 
Bernstein'^''  reported  French 
showed  a decrease  in  flow  stress  due 


are 


T2  - 1/2  sin  0 and  T3  = 1/2  S-^  cos  |-  (1  - 2v  + sin  |) 

where  are  the  principal  stresses,  with  and  S2  in  the  specimen  plane 
and  S3  ■ vCSi  + S2)  out  of  plane.  Figure  14  shows  Schijve’s  results. 

Since  subgrains  form  due  to  the  attainment  of  a critical  dislocation 
density,  and  dislocations  are  formed  by  shear,  the  stress  range  which  has 
been  calculated  and  is  shown  in  Figures  7 thru  10  is  probably  a shear  stress 
range.  The  asympotlc  solutions  which  have  been  derived  are  of  the 
wrong  special  form,  since  they  predict  a stress  singularity  at  the  crack  tip, 
but  their  angular  form  appears  to  be  correct.  Ahead  of  the  crack  tip, 

6*0”,  T3  > T2  gives  a finite  value  to  the  plastic  zone  in  that  region, 
which  T2  does  not.  Thus,  the  subgrain  forming  region  is  probably  due  to  the 
maximum  shear  stress  range  experienced  by  Che  material. 

It  is  not  possible  to  correlate  the  subgrain  size  distributions  as 
measured  and  graphed  in  Figure  1 with  those  of  other  investigators  because 
no  other  similar  data  have  been  found  covering  the  whole  plastic  zone. 
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There  have  been  several  transmission  electron  microscopy  investigations 

into  subgrain  formation  near  fatigue  crack  tips  in  iron  and  low  carbon 

steel which  do  indicate  that  subgrains  of  a minimum  size  of 

0.8  to  1.5  ym  are  found  adjacent  to  the  crack  plane  for  values  of  AK  * 15  MN/m3/2. 

At  higher  AK  (~  35  MN/m3/2)  smaller  subgrains  (~  0.5  ym)  may  form.^^l) 

The  special  form  of  the  crack  tip  stress  distribution  shown  by  Rice 
and  Johnson ^23)  their  Figure  4(a)  is  (as  pointed  out  by  Schwalbe (24)) 


O/Oy 


0.3 

0.1  + 7 


where  r » r/(AK/a)^. 


(26) 


This  equation  has  almost  the  same  form  as  equation  (19)  and  it  does  not  predict 
a singularity  at  the  crack  tip.  Rice  and  Johnson  determined  the  data  to 
which  equation  (26)  was  fitted  by  finite  element  analysis  of  a monotonlcally 
loaded  crack  for  a non-work  hardening  material;  the  similarity  between  their 
calculation  and  equation  (19) , which  was  -determined  experimentally  for  a 
cyclically-hardening  material,  is  thus  remarkable. 


(25) 

As  commented  on  by  Rice,  a small  volume  of  material  near  the  tip  ; 

of  an  unloaded  fatigue  crack  tip  experiences  compressive  stresses.  This 
point  was  further  amplified  by  Hahn,  et  al  (26)  their  etch  studies 

of  fatigue  crack  propagation  in  Fe  - 3 Si.  The  schematic  of  the  cyclic  i 

history  of  the  material  ahead  of  a fatigue  crack  tip  tip  (their  Figure  14)  { 

drawn  in  terms  of  strain  range  fits  well  the  result  expressed  by  equation 

(21) 


2 

A + Br  ' 


The  corresponding  expected  stress  response  of  the  material  is  shown  schematically 
in  Figure  15  for  both  the  dry  nitrogen  and  the  wet  air  environments.  For 
most  of  the  cycles  which  an  element  of  material  experiences  as  the  crack  approaches, 
the  stress  is  positive,  but  when  the  crack  tip  is  very  close,  the  large  strain 
amplitude  experienced  by  the  material  requires  compressive  stresses  to 
obtain  the  crack  closure  observed.  Since  the  strains  at  the  crack  tip  for 
the  wet  air  environment  are  less  than  those  for  the  dry  nitrogen  environment, 
lower  compressive  stresses  are  required.  The  existence  of  crack  tip 
compressive  stress  has  been  demonstrated  experimentally  by  Aira  and  Tanaka(22) 
who  found  an  Increase  in  amplitude  with  increasing  AK.  Since  their  measure- 
ments were  made  in  "laboratory  air",  these  observations  are  consistent  with 
our  finding  that  the  stress  range  Increases  with  Increasing  AK  for  the  water 
vapor  environment  (for  AK  up  to  about  13  MN/m^/2) . 

There  is  no  direct  evidence  for  the  assumption  that  a wet  air  environment 
lowers  the  compressive  stresses  in  the  unloaded  state,  as  cowared  to  the 
dry  environment,  but  there  is  indirect  evidence.  Jono,  et  al(2®)  have 
directly  observed,  at  high  resolution,  cyclically  loaded  crack  tips  in  grains 
of  Fe  - 3 Si  oriented  to  shear  only  In-plane  which  were  propagated  in 
vacuum  and  in  "laboratory  air"  environments.  In  the  vacuum  environment,  the 
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Figure  15.  Schematic  Representation  of  the  Stress  Range 
Experienced  By  Material  Ahead  of  the  Crack 
Tip  Showing  the  Effect  of  Environment. 
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crack  Is  found  to  open  first  behind  the  crack  tip,  with  the  crack  opening 
progressing  towards  the  tip  as  the  load  Is  Increased,  whereas  after  being 
cycled  In  the  laboratory  air  environment,  the  crack  begins  opening  at  the 
crack  tip.  The  authors  attribute  this  behavior  to  crack  "reweldlng"  In  the 
vacuum  environment,  but  our  Interpretation  of  this  finding  Is  that  the  crack 
tip  compressive  stress  magnitude  Is  greater  In  the  vacuum.  The  direct 
observations  by  Jono,  et  al  also  support  our  finding  of  a decreased  strain 
amplitude  at  the  crack  tip  due  to  the  presence  of  water  vapor  while  cycling. 
Careful  e^Mnatlon  of  their  crack  tip  photographs  by  a stereoviewing 
technlque^^’^  which  allows  visualization  of  crack  tip  displacements.  Indicates 
that  more  pronounced  shear  bands  with  higher  strain  magnitude  exist  for  the 
vacuum  environment  than  for  the  wet  air  environment.  The  crack  tip  Is  rounder, 
and  the  crack  tip  opening  displacement  (CTOD)  for  a value  of  K^gif  is  greater. 
Figure  16,  for  the  wet  air  environment.  The  change  In  hysteresis  of  CTOD 
with  change  of  environment  Implies  that  the  near-crack-tip  material  also 
undergoes  less  hysteresis  In  the  water  vapor  environment,  which  Is  In 
agreement  with  the  calculated  energy  for  crack  propagation. 

Several  Investigators  have  found  that  water  vapor  in  the  environment 
causes  a change  in  the  mode  of  crack  tip  opening.  Schijve^^®)  has  observed 
a change  In  mode  from  a predominantly  "shear"  to  a "tensile"  mode  on  going 
from  a vacuum  to  humid  air  environment  In  2024  and  7075  aluminum  alloys. 
Davidson  and  Lankford  previously  reported  a similar  result  for  Fe  - .05C 
steel.  These  observations  generally  agree  with  a decreased  size  of  the 
plastic  zone  for  the  crack  propagated  in  wet  air(^),  and  with  the  present 
finding  of  decreased  crack  tip  plastic  strain  and  stress  range  for  these 
conditions . 


The  Mechanism  of  Crack  Propagation 

From  the  data  so  far  gathered  In  this  Investigation,  coupled  with  the 
findings  of  other  Investigators,  the  effects  of  a water  vapor  environment 
on  fatigue  crack  propagation  In  low  carbon  steels  at  low  AK  are  as  follows: 


1)  The  mode  of  crack  propagation  is  changed  from  a mix  of  mode  II 
(shear)  + mode  I (tensile) , where  mode  II  predominates,  to  a mix 
of  cracking  Involving  more  mode  I opening.  The  shape  of  the 
crack  tip  probably  changes  also. 

2)  The  Intensity  of  the  shear  band  strain  Is  decreased  by  there 
being  more  bands  operative  having  lower  strain. 

3)  The  size  of  the  subcell  forming  region  Is  decreased,  and  the 
subcell  sizes  are  Increased  at  the  crack  plane. 


4)  The  magnitude  of  the  shear  stress  range  at  the  crack  tip  Is 
diminished.  The  proportionality  of  stress  range  to  r~^  Is 
approached  more  slowly. 


5) 


Stress  range  experienced  by  the  material  has  the  form 


Ao  - 2.26  o (A  + Br)'*^^^ 
o 
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Figure  16.  Crack  Tip  Opening  Displacement  As  a Function  of 
Stress  Intensity  Factor  For  a Crack  Propagated 
in  Humid  Air  va  Vacuum.  From  Photographs  in 
Jono,  et  al  Arrows  Denote  Loading  Path. 


and  strain  range 


Aep  “ 2/(A  + Br) 


where  values  of  A and  B are  given  in  Figure  1 for  both 
humid  air  and  dry  nitrogen  environments. 

6)  The  compressive  stresses  at  the  crack  tip,  for  the  unloaded 
condition,  are  probably  reduced. 

7)  The  energy  expended  in  producing  a unit  area  of  crack  is  strongly 
reduced  at  low  AK,  increasing  to  equal  that  for  the  dry  nitrogen 
environment  as  AK  increases. 

8)  The  presence  of  hydrogen  at  the  crack  tip  probably  a)  reduces 
the  Fe  - Fe  bond  strength,  and  b)  increases  the  cyclic  and 
monotonic  yield  strengths  of  this  low  carbon  steel. 

This  information  is  expected  (a)  to  be  useful  from  the  viewpoint  of 
mathematically  modeling  fatigue  crack  propagation,  and  seems  (b)  to  fairly 
well  delineate  the  mechanism  of  crack  propagation  and  its  alteration  by 
environment,  ^foreover,  it  is  possible  to  summarize  the  foregoing  concepts  into 
a fairly  cohesive,  albeit  somewhat  speculative  model,  as  follows;  although  there 
may  exist  alternate  explanations,  there  exists  very  little  current  information 
which  is  not  supportive  of  the  basic  ideas  developed  below. 

Fatigue  crack  propagation  (no  environmental  effect) 

At  low  values  of  AK,  fatigue  crack  propagation  is  exemplified  by 
mixed  mode  crack  opening  (mode  I + mode  II) . The  mode  I (tensile)  opening 
drives  the  crack  forward  while  the  mode  II  (shear)  displacement  dissipates 
energy  in  the  lattice  through  shear.  As  AK  increases,  the  amount  of  mode  I 
Increases,  driving  the  crack  even  faster.  Schematically,  this  may  be 
represented  as  in  Figure  17,  which  is  a modified  way  of  viewing  a similar 
concept  originating  with  Tomkins  and  Wareing.^^^i 

The  material  at  the  crack  tip  responds  to  the  mode  of  crack  opening 
by  shearing  in  one  or  more  shear  bands . The  number  of  shear  bands  and  the 
angle  of  the  shear  band  to  the  loading  axis  depend  on  a number  of  factors 
including  the  angle  of  the  cracking  plane  to  the  loading  axis  and  the 
magnitude  of  AK.  As  AK  decreases,  the  path  of  the  crack  becomes  much  more 
irregular,  so  as  the  threshold  is  approached,  the  energy  required  for  producing 
a unit  area  of  crack  (projected  onto  a plane  perpendicular  to  the  loading 
axis)  increases. 

The  environmental  effect 

The  major  effect  of  environment  is  to  alter  the  mix  of  opening 
modes,  causing  there  to  be  an  increase  in  the  mode  I (tensile)  magnitude  for 
a value  of  AK.  Since  the  shear  sustainable  by  the  material  is  lowered,  the 
number  of  shear  bands  is  Increased,  with  a decrease  in  intensity  of  strain 
in  each.  The  shearing  bands  also  have  a decreased  length  owing  to  their 
number.  CTOD  for  a fixed  strain  at  the  crack  tip  is  Increased  for  a value 
of  AK  because  of  the  Increased  mode  I opening  and  accompanying  greater 
number  of  less  intense  shear  bands. 
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Hydrogen  in  the  lattice,  which  may  come  from  the  catalyzed  dissociation  of 
H2O  in  the  presence  of  Fe  by  the  mechanism  described  by  Dwyer,  Simmons  and 
Wei(32)  probably  is  carried  into  the  lattice  by  a combination  of  diffusion^^S) 
and  dislocation  sweep  in^^^^.  The  effect  of  this  lattice  hydrogen,  then, 
is  probably  a combination  of  (a)  lowering  the  energy  of  the  Fe  - Fe  bond, 
making  the  rupturing  process  easier,  thus  requiring  less  mechanical  input  of 
energy,  and  (b)  Increasing  the  resistance  to  plastic  flow  in  the  material. 

This  latter  effect  concentrates  the  effect  of  strain;  thereby  reducing  the 
volume  of  plastically  deformed  material  and  preventing  some  of  the  strain 
at  the  crack  tip  in  individual  shear  bands.  Thus,  the  energy  required  for 
cracking  is  reduced.  The  reduction  of  strain  magnitude  in  each  shear  band 
decreases  the  requirement  for  reversed  shear  on  the  unloading  cycle,  which 
lowers  the  compressive  stress  at  the  crack  tip  in  the  unloaded  condition. 

This  implies  that  the  crack  tip  begins  to  open  at  a reduced  load  on  the 
tensile  portion  of  the  cycle,  so  that  a greater  CTOD  and  crack  growth  per 
cycle  would  be  expected. 
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